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■ Abstract. Critical points of a scalar quantitiy are either extremal points or saddle 

points. The character of the critical points is determined by the sign distribution of the 
eigenvalues of the Hessian matrix. For a two-dimensional homogeneous and isotropic 
£j ■ random function topological arguments are sufficient to show that all possible sign 

<4_H I combinations are equidistributed or with other words, the density of the saddle points 

c/2 ' and extrema agree. This argument breaks down in three dimensions. All ratios of the 

. ^ , densities of saddle points and extrema larger than one are possible. For a homogeneous 

Gaussian random field one finds no longer an equidistribution of signs, saddle points 
(— i ■ are slightly more frequent. 

(N ■ PACS numbers: 02.10.Yn, 02.40.Re, 47.54.+r, 47.27.Gs 

> ' 

t> : 

' 1- Introduction 

p ; 

1^! . Often turbulent phenomena are accompanied by random fluctuations of some scalar 
quantity. One may think of the density, temperature or optical refraction index for 
compressible flows, concentrations in the case of mixtures or reaction or combustion 
£> . rates in the case of chemically active substances. Important properties of these fields 

•i-H , 

are their critical points, i.e. those points where their gradient vanishes. They give 
rise to a subdivision of these fields into "dissipation elements" pQ very similar to the 
partition of the velocity field into "eddies" generated from the stagnation points [21 E]- 
One should also mention [I] where the routes to the generation of critical points out of 
simpler flows are studied. In addition to their importance for the characterization of 
the field, they are often also of direct physical significance. Let us discuss the situation 
in more detail in case of a temperature field. There the heat flux is directed towards 
lower temperatures, i.e. the heat flux lines end in a temperature minimum. Each 
temperature minimum is surrounded by a cell consisting of points with heat flux directed 
towards the minimum. It does of course not mean that the actual heat transport in 
the turbulent flow is appropriately described by this cell structure as the cells change 
from time to time and there usually is a significant contribution from convection to 
the heat transport. The cell structure is more suited to visualize the temperature field. 
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As the cells are associated with temperature minima the distribution of the minima is 
of special interest. The distribution of the minima depends of course on the random 
properties of the turbulence and is therefore difficult to obtain. Some hints can however 
be obtained if one assumes Gaussian random fields, an assumption which is often made 
in turbulence studies. The properties of the critical points depend on the Hessian 
matrix of the scalar field. The Hessian of isotropic homogeneous random functions 
has been studied for two-dimensional fields by Longuet-Higgins [5] in a paper on light 
reflection at a random surface and for three-dimensional fields by Halperin and Lax 
[6] in a paper on impurities in semiconductors. In these early works the method of 
Rice [7] to represent a random function as a superposition of trigonometric functions is 
used. This method assumes implicitly homogeneity of the random function, its role was 
therefore not considered explicitly in these works. More recently the Hessians of large 
dimensions were studied in the context of random matrix theory [8] by Fyodorov [91 [TO]. 
The relation of the Hessians to the Gaussian orthogonal ensemble (GOE) was clearly 
formulated in this work and it was shown that one can introduce a normally distributed 
auxiliary variable t such that the ensemble of Hessians of homogeneous isotropic random 
functions is obtained from GOE by averaging over this variable. Through this approach 
it is possible to apply results from GOE studies at the cost of one extra average (i.e. 
integration). Fyodorov shows how this can effectively be done. 

The topic which is of high interest to us is the distribution of extrema in the 
critical points. The type of the critical point is determined by the sign distribution of 
the eigenvalues of the Hessian. Usually the number of negative eigenvalues is called 
the index of the matrix. Methods of functional integration have been used by Bray 
and Dean [TT] to determine the asymptotic distribution of the index and Fyodorov , 
Sommers and Williams [12] have determined the asymptotic value of the minima i.e. 
the critical points of index zero. Comparing their eqs. (31) and (23) one finds that the 
probability of a critical point to be a minimum is asymptotically given by exp(-N) for 
N-dimensional homogeneous random functions. 

The situation is completelely different for two-dimensional random functions. There 
it is well known that it follows from topological reasons or from the vanishing of the 
average of the determinant of the Hessian that the densities of the extrema and of the 
saddle points agree. This is e.g. the basis of the Poincare- index and of the Brouwer- 
degree.[13j [T3] For random functions with even probability distribution, e.g. in the 
Gaussian case one has one half of the critical points with index 1 and a quarter with index 
zero or two. The four different sign combinations of the eigenvalues are equidistributed. 
On the other hand the result of Fyodorov et al shows that the fraction of critical points 
which are minima, namely exp(-N), is for large matrices well below the fraction of 2~ N 
to be expected for an equidistribution of signs. 

Our interest is the three-dimensional case. There the strong topological restrictions 
from the two-dimensional case no longer apply. From Morse theory (15] one can 
conclude that the density of the extrema cannot be smaller than the number of saddle 
points, but equality, as in two dimensions is no longer true. In an appendix, we 
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describe functions having various fractions of minima. So we consider Gaussian random 
functions. Although we could use results of Fyodorov for the critical point density, the 
situation is more complicated for the density of minima. So we stay with the ensemble 
of Hessians and use the methods developped in random matrix theory as far as they 
apply and make use of the additional symmetry which follows from the homogeneity 
of the random functions. We think that this symmetry can best be expressed in terms 
of the wedge product of the exterior algebra. One finds that the average of all wedge 
products of column vectors of the Hessian vanish. For the three-dimensional case one 
can then evaluate the remaining integrals by fairly elementary means and obtains for 
the fraction of minima under all critical points the value 1/4 — 3V/6/58 = 0.1233 slightly 
less than 1/8 to be expected for equidistribution of signs. 



2. Basic relations 



To obtain the density of the critical points of a scalar random function C(x), i.e. points 
in TV-dimensional cc-space where the gradient g = VC vanishes one observes that each 
critical point gives a contribution 1/| det(H(a;))| to the integral over 5(g), where 5(g) 
denotes the Dirac 5-function. Here H(x) denotes the Jacobian matrix of g, i.e. the 
Hessian matrix of the second derivatives of C and det(H) its determinant. Then the 
density of the critical points be obtained as 



ncr =<|det(H)|%)> with 9=[g^ )' H ={g^)- W 

where < a > denotes the average of a. This is the so-called Kac-Rice formula. Similarly 
the density n m i of the minima of C can be obtained from 

rVi=<Xmi(H) | det(H) |,%)>. (2) 

Here x m j(H) denotes the characteristic function of the minima which is one if all 
eigenvalues of H are positive and otherwise zero. 

From the homogeneity condition one obtains with integration by parts 

< HiiH kk >=< H? k > . (3) 
This is in marked contrast random matrix theory where matrices K from the GOE fufill 

< KaK kk >goe= 25 ik < K} k >goe ■ (4) 
The isotropy of C(x) lead with ([3]) to the probability distribution function (pdf) [9] 

P = cexp ( - ag 2 - (3( Tr H 2 - ^L^(Tr H) 2 )) (5) 

1 n 1 
with a = 5 — , p 



2<g(>" 2<H( 2 > 
where Tr denotes the trace of the matrix and c is a normalization constant. Obviously all 
off-diagonal elements are independant and have the same normal distribution, exactly 
as for GOE. Differences occur only in the diagonal elements. Therefore Fyodorov 



Critical point density 



4 



introduced an additional normally distributed variable t (averages over t are denoted by 
a subscript NO) and wrote with K from GOE 

H ik = K ik + t5 ik . (6) 

To determine the mean square of t one averages over GOE and NO to obtain for i ^ k 

< Hf k >=< HaH kk >=< (Ku + t){K kk + t) >goe,no=< t 2 >no (7) 
Alternatively it is possible to assume t = \ Tr K and one obtains for i ^ k 

< H u H kk >=< (Ku + X TTK)(K kk + X TrK) > GOE = 

2 X < K\ > GOE +X 2 N < K 2 U >goe=< H ik > 2 =< K\ k > GOE 

and therefore 

2N X 2 + 4 X - 1 = 0. 

This representation is useful for Monte Carlo calculations. 

Now we want to show that the homogeneity condition ([3]) implies that the average 
value of the determinant of H vanishes. Much more is true, namely the average of all 
coefficients of the characteristic polynomial of H vanish (with the obvious exception 
of the first one which is equal to one). This is in marked contrast to the GOE where 
already the average of the determinant is different from zero. This property of H can 
be conveniently shown using the exterior product of the exterior algebra [16]. This 
is a product h x A h 2 defined for p = 1,2, ... ,N on the vector space of p-tuples 
v ± = h x A h 2 A . . . A h p of vectors from some basic vector space of vectors hi which 
fulfills the associative and distributive laws and furthermore an anticommutative law 

h ± A h 2 = —h 2 A h ± . 

As basic vector space we take an N- dimensional vector space with basic vectors 
e ± .e 2 , ■ ■ ■ , c at and consider the row vectors of H, i.e. hi = Hne x + H i2 e 2 + . . , + HiNeN- 
Then we can write the homogeneity condition ([3]) as 

<hiAh k > = < (H u ei + H ik e k ) A (H ki ei + H kk e k ) > 
= < HaH kk — H ik H ki > e.i A e k = 0, 

i.e. the wedge product of two arbitrary row vectors of H has vanishing mean. As the 
average of an arbitrary number p of Gaussian random variables vanishes for odd p and 
is for even p just the sum of all partitions of the product into pairs and taking their 
averages one finds that all these averages vanish too, 

< h i% A /ij 2 A ... A h ip >= p = 1, 2, . . . , N. (8) 

These equations show the high degree of symmetry of this Gaussian matrix ensemble, all 
antisymmetric products have vanishing mean. One special case, namely p = N, states 
the vanishing of the average determinant, 

< det(H) > = 0. (9) 
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From (jSJ) one can obtain a much more general equation. With (jSJ) one has also 

< (h^ - dij A - a ia ) A ... A (h ijv - a ijv ) > = (-1)^0^ A A . . . A a iN 

with arbitrary constant vectors a p . This can also be written as 

< det(H - A) > = det(-A) = (-1)* det(A). (10) 

with an arbitrary constant matrix A. For A proportional to the identity matrix ffTOl) 
shows, that the average of all nontrivial cofficients of the characteristic polynomial of 
H vanishes. Related results for GOE can be found in [8j [17] where the moments of the 
characterstic polynomial have been determined. With ([6]) an alternative derivation is 
possible. 

An interesting application of (Q is possible in two dimensions. There the sign 
of det(H) is positive for extrema of C and negative for saddle points. Then we see 
from (|2j) that (Q implies equal densities of the extrema and of the saddle points in two 
dimensions, this argument, however, cannot be carried over to higher dimensions, there 
one has equal densities of critical points of even and odd indices (Brouwer degree). 



3. Three-dimensional matrices 



From (GO) one notices that g and H are statistically independent. So the density of the 
critical points and of the maxima can be written as a product of two factors n = n^ 2 \ 
where the first factor depends only on the average values of g and the second of those 
of H. The first factor is easily determined 



n 



(i) 



J 5(g)exp(-ag 2 )d 3 g 



Jexp(-ag 2 ) d 3 g 
It remains to find the second factor. 

Let us now restrict ourselves to three dimensions, N 
is then given by 

P = cexp(-/3(TrH 2 - (TrH) 2 /5)) 



3. The pdf for the Hessian 



with (3 



(12) 



2 < H( 2 > 

We are interested in averages depending only on the eigenvalues of H. Then it is 
common [5] to introduce the eigenvalues Ei,E 2 , E% as new variables and to average over 
the remaining variables. As the ordering of the eigenvalues is arbitrary, we assume 

E 1 <E 2 < E 3 . (13) 

Then it is shown in j8] that the pdf in the eigenvalue space P e is given by 

1 



P e = cexp(-f3e 1 Ae)A(e) with A 



ik 



-(55, 



where A(e) = (E 3 — E 2 )(E 2 — Ei)(E 3 — E\) and e is a column vector of the eigenvalues. 
Now we introduce new variables 2/1,2/2, 2/3 through 
( yl\ , ( -2 Vb \ 



B 



2/2 
V j 



1 

V6 



V 



-y/3 V$ 

\/3 ) 



y 



(14) 
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Figure 1. The unit sphere in y-space. Bright is the physically significant region 
(|13p . Shown is also the partition into regions of only positive eigenvalues (topmost 
region), one negative eigenvalue {Ex < 0), two negative eigenvalues (E2 < 0) and all 
eigenvalues negative. 

One verifies easily that B T AB is the identity matrix. Therefore the pdf reads in the 
new variables 

P e = cexp(-(3y T y)A(By). 

With spherical polar coordinates 

yi = r cos ip sin yi = r sin ip sin 2/3 = r cos d 

one obtains A(e) = r 3 A 3 ($, cp) = r 3 sin 3 $ sin (3<p)/y/2. As we are interested in the 
average of homogeneous quantities $(e) = r n $ n (-$, <p) only, we can write 

< <3> > = J r n+5 exp(-/3r 2 )$ n ($, <p) sin(3^) sin 4 $ drdy?d$. 
Now we can perform the r-integration and obtain 

< $ > = ^tfn^ / ^ *P) Sin4 ^ ^ 

where T denotes the T-function. 

Now we have reduced the integral from one in the 6-dimensional space of the Hn, 
to an integral over a 2-dimensional sphere. One finds easily, that E\ = E2 corresponds 
to p = I and E2 = E3 to (p = 0. Then, because of the restriction ([TBI the integration 
is not over the whole sphere but only over the segment 

0<y?<7r/3, < ■d < 7T. (16) 

The other regions of the sphere correspond to a different ordering of the eigenvalues. As 
the enumeration of the eigenvalues is usually arbitrary, one can as well integrate over 
the whole sphere and divide the result by the number of different orderings, namely 
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six. The configuration is shown in figure 1 where we additionally show the cross section 
of the sphere with the planes E\ = 0, E 2 = 0, and E 3 = 0. As these planes, as well 
as Ei = Ei and E 3 = E 2 , all pass through the origin, the cross sections are all great 
circles. The segment (1131) is then divided into four spherical triangles. The topmost is 
bounded in addition to (fT6l) by E\ = 0, it corresponds therefore to eigenvalues which are 
all positive, the next triangle is bounded by E\ = and E 2 = 0, corresponds therefore 
to one negative eigenvalue, the next one to two negative eigenvalues and the last one 
to eigenvalues which are all negative. One also notices that the segment and the pdf is 
not changed if the signs of the eigenvalues and their numeration is reversed. The four 
spherical triangles are however interchanged. Therefore the first and the fourth and the 
second and the third lead often to contributions of the same size. 



4. Critical points 

As a first step, we determine the normalization constant c of the pdf. As usual, it can 
be determined from the average of one which is one. As one is a homogeneous function 
of degree n = one has from (fT5l) 



K ) 1 sin(3</?)dy? / sin 4 tfdtf = — ^U. (17) 



a/2/3 3 A) K ^' Mo 2 a/2/3 3 ' 

Now we determine the second factor nS 2 ' for the spatial density of the minima. As the 
determinant of H is a homogeneous function of third degree we obtain from (fT5j) 

n$ =< Xmi(H) det(H) >= J det(H) sin(3^) sin 4 dipcM (18) 

which is to be evaluated over the unit sphere in y-space where all eigenvalues are 
positivjjl, i.e. the topmost triangle in figure 1. With 

d*(H) = r'* 3 (^) = (^(5 -8^*) - ■ 
and the indefinite integral 

1 = J $ 3 (^)A 3 (^) sintfdtf = h%<p)+I 2 (#,<p), 
15 sin 3f /sin 5 $ 2sin 7, # N \ 



9 V 4 7 



T , . . a/3 sin 609 cos $ ( 8 4sin 2 $ sin 4, (9 sin 6 ^\ 

U$ m) = 1 L 1 1 1 (19) 

2V 6 V 105 105 35 42 / K J 

we obtain for the integral I m i over the topmost spherical triangle 

tt/3 tt/3 

// mi =r (i(& 1 (< P ),<p)-i(o,<p))d<p= [ (i(m<p),<p))*<p (20) 



I One might think that this can easily be done with a computer algebra program. In our experience 
these programs nead a lot of assistance which is indicated by the following intermediate results. They 
are nevertheless of enormous value in confirming the results. 
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as the integral at i? = vanishes obviously and where describes the great circle 

E 1 = 0, i.e. from flE 



2COSV9 

Smi7i = ,. , COSi7i 



a/5 + 4 cos 2 if ' \/5 + 4 cos 2 <£> 

Now the indefinite integrals IIi((p) and IIi((p) of Ji and J2 over E\ — can again be 
determined easily, and one obtains 

rr / ;\ v^cosfw) / 15 n 2/ . 356 cos 4 (09) \ 

= 5 ( - T + Scos 2 ^) — 

189^5 + 4 cos 2 (^) 4 25 

, ^""z ,(l-^- (21) 
9J5 + 4cos 2 (( / 9) 25 



816cos%) 3328cos 6 (y>) 2048 cos%) 

875 + 2625 + 2625 
With these indefinite integrals we get for II = II\ + II2 

//(-) = - — , 11(0) = -^, J/ mi = — (-v/3-V2) (22) 
v 3 ; 210' V ; 3780 ' 210^18 ' y ' 

/q\ _____ 

and then also n^l from (fT8l) where we have used the nomalization ([17]) 

"-=16^(1^"^) 
and these gives finally 

for the spatial density of the maxima. 

To determine the probability of the saddle points, we have to integrate over the 
spherical triangle between the great circles E\ — 0, E 2 — and E2 = E3. The first one 
is given by •& = $i(<£>) and let the second be given by 1? = $2 (<£>)• From (fT4"l) we obtain 

v 7 ^ 2cos(^ + tt/3) 

smi/o = =. cos Wo = - 



/5 + 4cos 2 (v9 + vr/3) ^5 + 4cos 2 (^ + tt/3) 

Then the integral over this triangle II sa can be written as 

I/ sa = / (J(#i, ¥>) - J(# 2 , V)) dy? = J/ mi - f (J(tf 2 , </?)) d^ (24) 



where we have taken into account, that the integrand is negative between E± = and 
E 2 = 0. Now p = <p + 2tt/3 gives tf 2 (<p) = $i(<p) and I(& 2 (<p), V 5 ) = and 
therefore 

r i{^ip)dip=r i(# 1 (<p),<p)d<p=- 1 i{$ 11 $)&{p=2ii{\). 

JO J-2tt/3 3 

We then obtain from ([21]) and fT22j) 

//. = — + 

210 V18 ' 

and for the ratio of the density of minima to that of critical points 

Una = I _ = 1233 

2II mi + 2II sa 4 58 
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5. Conclusion 

The distribution of critical points of a two-dimensional function is severely restricted 
by topological considerations. There the density of extremal points and saddle points 
agree, i.e. there is an equidistriburion of signs of the eigenvalues of the Hessian. For 
three-dimensional functions these topological restrictions no longer apply and this ratio 
depends on the type of the function. For a Gaussian homogeneous and isotropic random 
function sightly less than 1/8 of all critical points are minima. A comparison with results 
from computational fluid dynamics seems interesting. 
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Appendix: Example functions 

A simple function with many critical points is given by 

COS 1XX\ + COS 71X2 + COS ITX3. 

Its critical points are the lattice points (xi, x 2 , x 3 ) = (ni, n 2 , n 3 ) with integer ni, n 2 , n 3 . 
Maxima occur when the integers are all even, minima when they are all odd. All other 
critical points are saddle points. Obviously all sign combinations are equidistributed. 
Under eight critical points one finds one maximum one minimum and 6 saddle points. 
We now indicate how functions with arbitrary ratios not smaller than one of densities of 
saddle points to extrema can be found. Actually one can conclude from Theorem 3.4 in 
the Morse theory chapter of [15] that this restriction on the critical points is a necessary 
one. There are no functions having a density of extrema larger than the density of 
saddle points. 

The purpose of this appendix is to describe functions with various ratios of saddle 
point density to extrema density. As a first step, we describe that it is under certain 
conditions possible to combine two functions to a new function which has exactly the 
critical points of both combined. To be specific, let F(x) and $(cc) be two smooth 
functions with F(xq) = and $(oo) = and furthermore 

\VF(x) \ < M\x - x \ and |V$(x)| < 

then with some nonvanishing vector V the functions 

/ = a + V ■ (x - x ) + F(x) and 4> = V ■ x + $(x) 

are such that / has no critical points in a neighbourhood of Xq and <fi has only a finite 
number of critical points. Let a(r) be a smooth non-increasing function which is one 



\x\ 
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for r < 1 and vanishes for r > 2. Then one can define with some positive e 

9 = a(^) (. + + (1 - 

= a + y • (x - x ) + ( £ '<|>(!iZ£°) - F^)) + F(x). 

It is obvious that g has for sufficiently small e in |x — a; | < e just the critical points of 
and for |x — x | > 2e just those of /. With standard estimates one can verify, that 
for small enough e there are no critical points in the transition region e < \x — x \ < 2e. 

Now we describe functions having an arbitrary non-zero even number of saddle 
points and just two extrema. These functions can then be inserted in arbitrary numbers 
in some given smooth function and one can use them to generate functions having in any 
given domain an arbitrary rational ratio of saddle points to extrema above one. Then 
it is also possible to build functions with a similar density ratio. These functions are 
necessarily three-dimensional. We nevertheless begin with a function of two variables 
r, z, namely 

f(r, z) = z + a ( 1 + 6(r2 + (2 + 1) 2) " 1 + 6( r 2 + - l) 2 ) ) 

with positive a and b. Later on we will think of r and z as cylindrical coordinates in 
a three-dimensional space. The discussion of this function is elementary, so we just 
describe the results. For small values of a one has only a small perturbation of the 
linear function z without any critical points. For somewhat larger a there are four 
critical points, one maximum, one minimum and two saddle points, all situated on the 
z-axis r = 0. For even larger values of a, exactly for 

4a > (b + ^ 

the extrema remain on the z-axis while the saddle points move on the r-axis to r = ±r 
with some positive r s and this is the case which is of interest to us. In three-dimensional 
space with cylindrical coordinates this function has however a non-generic critical circle 
r = r s , so we would like to replace a by acosm?, d being the angle in cylindrical polar 
coordinates. This function is not even continuous at r = 0, therefore we modify a to 
some function a(r 2 ). It is easily checked that the number and type of singularities is not 
changed if a(r 2 ) is a monotonically decreasing positive function of r, which we assume. 
So, we have with 

, / q\ __ . 4fe(q (r 2 ) + ai(r 2 )cos(m9)) 

J{r,Z, )-Z + ( 1 + b ( r2 + ( z + 1 y^( 1 + b ( r2 + ( z _ 1 y^ 

a smooth function which has a maximum and a minimum, both on the z-axis and 
In saddle points in the z = plane, provided a (r 2 ) + a 1 (r 2 )cos(n , &) is for all d a 
monotonically decreasing positive function with 

4a (0)) > (b+l 



s 



b, 

and ai(r 2 )/r n remains bounded at r = 0. 
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